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ABSTRACT 


The  exact  solution  of  the  scattering  of  obliquely  incident 
plane  waves  by  an  elliptical  dielectric  cylinder  is  obtained.  It 
is  found  that  each  expansion  coefficient  of  the  scattered  or 
transmitted  wave  is  coupled  to  all  coefficients  of  the  series  ex¬ 
pansion  for  the  incident  wave  except  when  the  elliptical  cylinder 
degenerates  to  a  circular  one.  Both  polarizations  of  the  inci¬ 
dent  wave  are  considered:  one  with  the  incident  electric  vector 
in  the  axial  direction,  and  the  other  with  the  incident  magnetic 
vector  in  the  axial  direction.  It  is  noted  that  in  the  general  case 
of  oblique  incidence  the  scattered  field  contains  a  significant 
c ros s -polarized  component  which  vanishes  at  normal  incidence. 
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I.  INTRODUCTION 


The  problems  of  scattering  of  waves  by  a  circular  dielectric 
cylinder  have  been  considered  by  many  authors  (Refs.  1-3).  Most 
recently,  the  exact  solution  of  the  problem  of  diffraction  of  normally 
incident  plane  waves  by  a  dielectric  elliptical  cylinder  was  obtained 
(Ref.  4).  However,  the  general  case  for  oblique  incidence  has  not  been 
considered.  It  is  the  purpose  of  this  paper  to  present  a  complete  solution 
for  this  general  problem.  It  is  hoped  that  the  results  will  be  applicable 
to  the  problem  of  scattering  of  light  by  noncircular  fibers. 


II.  FORMULATION  OF  THE  PROBLEM 

To  analyze  this  problem,  the  elliptical  cyclinder  coordinates 
(£,77,  z),  as  shown  in  Fig.  1,  are  introduced.  In  terms  of  the  rectangular 
coordinates  (x,y,  z),  the  elliptical  cylinder  coordinates  are  defined  by  the 
following  relations: 

x  =  q  cosh  £  cos  Tj 
y  =  q  sinh  £  sin  7) 

(D 

z  =  z 

(0  <£  <  co,  0  <  rj  <  2tt) 


-1- 


Fig.  1  Plane  Wave  Incident  On  An  Elliptical  Dielectric 
Cylinder.  Arrow  Indicates  The  Direction  Of 
Incident  Wave. 


where  q  is  the  semifocal  length  of  the  ellipse.  The  contour  surfaces 
of  constant  £  ,  are  confocal  elliptic  cylinders,  and  those  of  constant  rj 
are  confocal  hyperbolic  cylinders.  One  of  the  confocal  elliptic  cylinders 
with  £  =  £0  assumed  to  coincide  with  the  boundary  of  the  solid  dielec¬ 
tric  cylinder,  and  z-axis  coincides  with  its  longitudinal  axis.  It  will  be 
assumed  that  this  dielectric  cylinder,  which  has  a  permittivity  a 
permeability  |±jY  and  a  conductivity  of  zero,  is  embedded  in  a  homogeneous 
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perfect  dielectric  medium  ( €0,  p0).  A  possible  solution  of  the  wave 
equation  is  then  R(£)  6(17)  e^Z  e“lwt  where  R  and  0  satisfy  the  differ¬ 
ential  equations 

~  "  -  (c  -  Zy^  cosh  Z£  )  R  =  0  (Z) 

d£2 

+  (c  -  2y2  cos  27;)  0  =  0  (3) 

drf 

in  which  c  is  the  separation  constant  and  \2  =  (k2  -  P2 )  q2/4,  k  being  the 
wave  number. 

The  periodic  solutions  of  the  equation  in  7]  are  of  two  types:  even 
about  ^  -  0,  and  odd  about  =  0  (Ref.  5).  They  are  possible  only  for 
certain  characteristic  values  of  c.  The  even  and  odd  functions  are  re¬ 
spectively  denoted  by  cen(rjf  and  sen(77#  y2)  with  the  sequence  in  n 
according  to  increasing  values  of  c .  It  should  be  noted  that  these  functions 
are  orthogonal  functions.  The  solutions  of  ( Z)  corresponding  to  the  even 
function  cen(7),  y2)  having  the  same  characteristic  values  of  c  are 
Cen  (6  Y2)  and  Feyn(£,  y2),  and  those  corresponding  to  the  odd  function 
sen(T7*  V2)  are  Sen(£,  y2)  and  Geyn(£,  y2)  (Ref.  5). 

The  proper  choice  of  these  functions  to  represent  the  electro¬ 
magnetic  fields  depends  upon  the  boundary  conditions.  For  the  region 
within  the  dielectric  cylinder,  all  field  components  must  be  finite.  All 
field  components  for  the  scattered  wave*  must  satisfy  the  Sommerfeld’s 

‘  ‘  '■ .  '  ;  ®  v ;  ■/*'" 
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radiation  condition  at  infinity.  Consequently,  the  appropriate  solutions 
of  the  wave  equation  for  the  region  inside  the  dielectric  cylinder  are 


Y  j  )  cen(77*  yf) 

lSen(C»  vf)  sen(i7,  Vi> 


i^z 


and  those  for  the  scattered  wave  are 


'"MeJ1I)>  (2)(£,  Yq)  cen^V’  Yq)  ^ 


>  e 


ipz 


where 


UeL1),(2)(e.  Y2o)  sen(^,  Yq)  J 


Me{l),{2)^f  y2)  8  Cen(^j  Yo)±iFeyn(^,  y*) 


(  1  )  ,  (2)  O  O  O 

Nen  (£>  Yo)  =  Sen£.  yJ)  ±  i  Geyn(£,  YJ) 


and  are  given  in  the  next  section. 


(4) 


(5) 


(6) 

(7) 


III.  SCATTERING  OF  AN  OBLIQUELY  INCIDENT  PLANE  WAVE 

Two  types  of  incident  waves  are  possible.  The  one,  called  an  E- 
wave  is  defined  by  Hz  =  0,  and  the  other  called  an  H-wave  is  defined  by 


Ez  =  0. 


.  .  ‘.S  i-yl 


It  can  be  shown  that  a  plane  wave  with  its  direction  of  propagation 

defined  by  the  angles  <j£>  and  0  (see  Fig.  1)  is  given  by 

ik0(x  cos^>  sin  0  +  y  sin  sin  0  +  z  cos  0) 
e 

ikQ  |^q  sin  0  (cosh  ^  cos^  cos  +  sinh  £  sin  T)  sin  )  +  z  cos0  j 


oo  r  i 

“  p?  ^e2n  £  *  V0  )  ce2rSV>  VG)  ce2n  VQ) 

n=  o  L 


7777  Se2n+2  <£•  Vo)  se2n+2  Yo)  se2n+2  <<£>  Y2) 

Zn+  Z 


+  P2n+1  C62n+1  (^'  Yo)  Ce2n+1  <77'  Yo)  ce2n+l  Y*) 


(8) 


+  7777  SC2n+1  (^’  Y°}  Se^+1  (77’  se2n+ 1  (<£>  Yq)] 


2\1  e^°Z  COS  ^ 


where  =  k^q^  sin^0/4,  k^  =  2tt/X0,  and  XQ  is  the  free-space  wave¬ 
length.  p2ni  P2n+1'  s 2n+ 2*  ands2n+l  are  joining  factors  (Ref.  5). 

In  order  to  simplify  the  notations  for  the  Mathieu  and  modified 
Mathieu  functions  without  any  ambiguities,  the  following  abbreviations 
are  used: 


Cen<£'  Yo  >  =  Cen(£) 
Sen  (£»  Yo)  =  Sen(£) 
Cen({,  yj)  =  Ce*(^) 

Sen(e,  Y2)=Se*(£) 


cen('7>  Yo>  =  cen^V) 
sen  i7]  >  yh  =  sen(7}) 

cen(r)>  y \ )  =  cen(v) 

sen(ij,  yf )  =  se*(i7) 
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o 


Men(1)’(2)  (£,  Yo)  =  Men(1),(2)  (£) 


(1).(2) 


(D,(2) 


Ne_ - <e.Yo>  =  Nen'-"'-'£> 


(9) 


with  Yj  =  (kj  “  ko  cos^  9)  q^/4  where  k^  = 


A.  E  -WAVE 

The  axial  components  of  an  incident  E-wave  are 

Ez  =  Eo  {  t^ie  rigkt-hand  side  of  Eq.  (8)|  (10) 

H1  =  0.  (11) 

It  is  interesting  to  note  that,  unlike  the  case  for  a  normally  incident  wave 
or  for  a  perfectly  conducting  cylinder  at  oblique  incidence,  the  boundary 
conditions  for  the  present  general  case  of  oblique  incidence  on  a  dielectric 
elliptical  cylinder  cannot  be  satisfied  if  the  z  component  of  the  scattered 
or  transmitted  magnetic  field  is  taken  to  be  zero  for  an  incident  E-wave, 
or  if  the  z  component  of  the  scattered  or  transmitted  electric  field  is  taken 
to  be  zero  for  an  incident  H-wave.  Hence,  referring  to  (4)  and  (5),  we  see 
that  the  axial  components  of  the  scattered  field  and  of  the  transmitted  field 
inside  the  dielectric  cylinder  must  be  of  the  form 
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Z'z  -  2E0  I  [fe  M4n<«>c*2n<^>  “an'*’ 


n-o-P2" 

B2n+2  ..  (') 


*2n+2  Ne2n+2(^  Se2n+2(,?)  8e2n+2(^>) 


CO 

HI  =  2E°n?0  Me2n^)ce2n(’7)ce2n(«/» 


+  "^f  Ne2l  +  2l°8e2n  +  2(7?)  “20+2  «*» 


+  1  ^TMe2"  +  l  Ce2n  +  Il)?)  ce2n  +  I  ^ 


+  i 


.  D2n+|  „  U) 


2nd-  I 


Ne2n+I  (C,se2n  +  l  (77>  Se2n-H  (^]  e'k°2 


cos  6 


and 


EJ  =  2Eo  X  [  -p!1-Ce2n(^)ce2n(^)  ce2n(<fr 

.  n=0  L'  2n 


+  Se2n  +  se*2n  +2 {7))  se2fH-2(</>) 

.  2n  +  2  .  • 


(12) 


(13) 


+ 1  Ict1  S*2"+'(JI  "L,™  “2„+,H  * 


ik0z  cos  6 


2n  + 


(14) 

Hz,  =  2EoZ  [  {°  Ce2n<7?)  Ce2n(<£> 

n  =  0  2n 

+  ^S‘*2n«<«'”2„  +  2W!ea.+  2^1 

+  i  Ce2n  +  Il?lce2r,  +  I(7'1  ce2-.+  l(^ 

2n  +  l 

+  ‘^77s4+l,^^+,l’'1“2n+,H*lk»ICM9 


(15) 


where  An,  Bn,  Cn,  Dn>  Fn,  Gn,  Pn,  and  are  arbitrary  unknown 

coefficients  that  can  be  determined  by  applying  the  boundary  conditions. 

*  *  *  ,  *  .  .  . 
p2n*  P2n+1'  s 2n+  2*  and  s2n+l  are  Joimng  factors.  All  transverse  fields 

can  be  derived  from  Maxwell’s  equations  with  the  knowledge  of  the  axial 

fields . 

.  The  boundary  conditions  require  the  continuity  of  the  tangential 
components  of  the  electric  and  magnetic  field  at  the  boundary  surface 
£  =t0;  i-e.. , 
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j^r.h.  s.  of  Eq.  (10)  with  C  =  CQ]  +  jjr.  h.  s.  of  Eq.  (12)  with  C  =  CQ] 

=  [r.h.  s.  of  Eq.  (14)  with  C  =  C0]  ,  (16) 

.  h.  s.  of  Eq.  (13)  with  £  =  £0j  =  Jr.  h.  s.  of  Eq.  (15)  with  £  =  CQ]. 


(17) 


^r~2  e  {ik°  cos0^[r-h-s-  of  Eq-(10)  with  C=  Ca] 


+  ikQ  cos  0  [r.  h.  s.  of  Eq.  (12)  with  £  =  CQ] 


-  o  .  h .  s .  of  Eq .  ( 1 3 )  with  C  =  C  Q  j 


(k^  -  cos ^  0) 


|ik0  cos  0  j-  ^r.  h.  s.  of  Eq.  (14)  with  C  =  £0] 


■iw^l  jg-  ^r-h.  s.  of  Eq.  (15)  with  C  =  j|. 

— ; —  {"^oaTT-  fr*h*s-  of  Eq.  (10)  with  £  =  £  1 

ih^0  L  u«o  L  °J 


(18) 


k**  sin  0 
o 


-HJ 


[r.h.s.  of  Eq. (12)  with  £  =  CQ] 


■iko  cos  0  —j  £  r .  h.  s .  of  Eq.  (13)  with  C  =  J} 


{  _iwei  j|~  [r- h- s*  of  Eq-  (14)  with  C  =  Co] 


(k2-k2  Cos20) 


-ikQ  cos  0 ^  £r.  h.  s  of  Eq.  (15)  with  £  =  Co  ]  > 


(19) 
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where  r.h.s.  means  the  right-hand  side.  It  is  noted  that  in  contrast 
with  the  spherical  or  circular  cylinder  case  the  angular  functions  in  the 
elliptical  cylinder  case  are  functions  not  only  of  the  angular  component 
but  also  of  the  characteristics  of  the  medium.  Consequently,  the  sum¬ 
mation  signs  and  the  angular  functions  in  the  above  equations  may  not  be 
omitted.  However,  it  will  be  shown  that  this  difficulty  may  be  overcome 
by  the  orthogonality  properties  of  Mathieu  functions.  Substituting  the 
expansions 

00, 

<**(77)  =  £  am  »n  (20a) 

n  -0 


00  . 


sem(7?)  =  £  Pm  n  *en  t77J 

n  =0 


d 

d77 


00 


cem(T?)  =  I  rm  n  sen 

n-0 


{r/) 


00, 


sem(77)  =  X  Xm  n  ce  (77) 
a77  n  =0 


(20b) 


(20c) 


(20d) 


into  Eqs .  ( 1  6)  through  ( 19), .  and  applying,  the  orthogonality  relations  of 
Mathieu  functions,  leads  to  the  following  expressions: 


p*  [ Cen(Co>  +AnMe^)(f0)l  cen(<£)  =.£  -dr  ^m^o1 
n  L  J  .  m=0  m 


«£) 


*m,n 


(21) 
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* 


Sen(£0)  +BnNen(Co)]  sen(<£)  =  X  T*Senfio'**eTn(4>  ]  Prn,  r\  ' 


(22) 


%  Mb^0)  cen(4>)  =  X  -p*  Cem  .(C0)  cem  (^}  am,n 

.m=o  m 


(23) 


Neln  (£0)s en(<jb)  =  X  JIT  Se^ (£0 ) sem(<£)  /3m,n  • 

s"  m=o  m 


(24) 


5~cosfl(l-  ft  Si2n2C)  X  J-[sem(Co)+BmNe(^o)]sem^)Xmin 
/^o  \  k  j  -  kQcos  a  /  tti=o 


^-MenU/(C0)cen  «£) 


f-g-l)  ( J° ff-V)  I' -ft)  cSMc.mtt>  •»,,  . 

\  /^o/  \  k,  -  k«  cos  a  /m=o  171 


(25) 


II  COS0  (l-  kg°-K^e)  £  -pS,[Ce^^o)+A^><fo)]  cem«£>7m,n 


-  xr  Ne^  ^  (£o)  8en  (*£) 
5pi  n 


(26) 
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"g~  Sen  (^0)  +  Bn  Nen  (£0)  sen(<£ 


/M 


£„  m." 


k§  sin20  «,  Gm  ;  j0 

kf-k*cos40  8m*  Sem  ^o)sem^)/3m,n  > 


(27) 


Pn 


[ce,-  (Co>  +  An  Men  (£o>]cen(<£) 


+ 


COS 


k|  sin20  \  ?•  Dm 
kf-k2cos2  0/  £0  s" 


Kq  sin2  0 

kf  -  k|.  cos2© 


(C0)cem(^> 


m,  n 


(28) 


(n  =  0,  I, 2,3  ,4  .  .  .  .  ). 

am,rt,£m,n»ym,n*and  Xm,n  are  given  in  the  Appendix.  The  prime  on  the 
summation  sign  means  that  when  n  is  odd,  the  above  series  are  summed 
over  all  odd  values  of  m,  and  when  n  is  even,  the  series  are  summed 
over  all  even  values  of  m.  The  primes  on  the  modified  Mathieu  functions 
denote  differentiation  with  respect  to  £0.  The  unknown  coefficients  An# 

^ n,  Gn>  Pn,  and  Qn  can  now  be  obtained  from  the  above 
equations.  Combining  Eqs.  (21 ),  (24)#  (26),  and  (28)  gives 


0-. 


e  ■.  e 
Z*  ^ms  mn  “**  ^rri  *  mn 
m~0 


Z  Fmumn  +  Qmvmn=wn 

m«o 


(  n  =0  }l  ,2,3  .  .  .  ) 


and  combining  Eqs.  (22),  (23),  (25)  and  (27)  gives 


o°  / 

X  ^msmn 


-f  p  t 
~  rm  1 


0  =0 
mn 


m=o 


* 


(29) 


(30) 


(31) 


o  o 

/.  ®mumn"FFmvmn  ' 
m=o 

(n  =  0,  I,  2,3  .  .  .  ) 


(32) 


where 


_  00  / 

e  _  /£o  C0SQ  ( |  -  x2)  f  Z  amr  7mn 
smr>  V/io  m  r% 


00/ 


cos#  U-*2)flmZ 
r*o 


^mrXmn 


tmn= 

tmn=  [ /it  x8hfn  "  hm]  amn 


13- 


#  Umn=.[fm‘45",tZ  "ft  'fk''' ' 


%  ,  *  Vi 

!»-V  ■. v  • 

lJ'  ■>  ■/* 

bn  ■  »■;.€■  I 


.  0  ....  r  ®n  I  /  1  n 
.^rnn  =  I  9m  g  9m  J  "mn 


mn 


CD  / 


L  =  a/t^cosSi,-x2)^I  /3mrXrn 


6  o 


r  =o 


o 

vmn" 


oo  / 

cos#(l-x2)  fm  X 
r  =  o 


frmr 


Xrn 


/ 


e  bn'  , 

Wn  =  9n“b^  "9" 


(33) 


The  following  abbreviations  have  been  used: 


hn  =^-  Cen(Co)  cen(<p), 


hn  =  ~  Cen((0 )  cen(<p), 
Pn 


in  =  ~  Sen(Co)  sen 
sn 


Sn - ^en(£o)  sen  (</>)> 


=  -L  uBU) 


= 

n  Pn 


Me},  '(Co>  cen^>- 


=  -L  mJD’ 


an  =  Men  <*o>  cen(^)- 


1  (1) 

bn  =  —  Nen  (CJ  sen (</>). 


-  1  (1)' 

bn=—  Nen  (Co)  sen(^>), 
sn 
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The  coefficients  Fm  and  Qm,  and  .Gj^  and  Pm  can  be  obtained  readily 

from  Eqs.  (29)  and  (30),  and  Eqs.(31)  and  (32)  respectively.  Knowing 
Fm,  Qm*  Qm»  anc*  coefficients  An,  Bn,  Cn,  and  Dn  can  be 

found,  respectively,  from  Eqs.  (21)  through  (24). 

It  is  interesting  to  note  that  the  roots  of  the  determinant  of  Eqs. 
(29)  and  (30)  provide  the  propagation  constants  of  a  set  of  surface  waves 
along  an  elliptical  dielectric  cylinder  (Ref.  6).  Due  to  the  asymmetry  of 
the  elliptic  cylinder,  it  is  possible  to  have  two  orientations  for  the  field 
configurations.  The  propagation  constants  of  the  other  set  of  surface 
waves  are  obtained  from  the  roots  of  the  determinant  of  Eqs.  (31)  and  (32). 
This  surface  wave  problem  has  a  direct  application  to  the  problem  of 
guiding  light  waves  along  a  flat  fiber  i.  e.  ,  the  fiber  optics  problem 
(Ref.  6)  . 

The  transmitted  fields  in  the  dielectric  cylinder  and  the  scattered 
fields  in  free  space  due  to  an  incident  E-wave  are  now  completely 


determined. 


B.  H-WAVE 


The  corresponding  result  for  the  case  of  an  incident  H-wave  is 
of  the  same  form  as  above  if  E  is  replaced  by  H,  H  is  replaced  by  -E, 
€  is  replaced  by  j and  p,  is  replaced  by  e  ,  throughout. 


IV.  CONCLUSIONS 


The  exact  solution  of  the  problem  of  the  scattering  of  an  obliquely 
incident  plane  wave  by  an  elliptical  dielectric  cylinder,  or  by  a  dielectric 
ribbon  is  obtained.  It  is  noted  that  unlike  the  case  for  a  circular  dielectric 
cylinder,  each  expansion  coefficient  of  the  scattered  or  transmitted  wave 
for  the  elliptical  dielectric  cyclinder  is  coupled  to  all  coefficients  of  the 
series  expansion  for  the  incident  wave.  Hence,  the  results  are  much  more 
involved.  Numerical  investigation  for  the  case  of  normally  incident  plane 
waves  shows  that  the  infinite  determinants  representing  these  expansion 
coefficients  for  the  scattered  or  transmitted  wave  converge  quite  rapidly 
for  small  values  of  kQq/Z;  only  the  first  few  terms  in  the  determinants 
are  needed  as  long  as  k0q/2  is  less  than  5  (Ref.  4).  This  range  of  kQq/2 
covers,  however,  most  of  the  range  not  covered  by  the  usual  approximate 
diffraction  theory.  The  results  presented  here  are  particularly  useful  in 
studying  the  scattering  of  light  by  thin  fiber  ribbons. 
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.  APPENDIX 


FORMULAS  FOR  atn>  f*m<  n,  ym>  h>  AND  Xm,  n 


It  can  readily  be  shown  from  the  theory  of  Mathieu  functions  that 


27T  -2  rr 

otm,n  =  f  ce*  (17)  cen  (t?>W/  ce*  (77)^17 
...  Jo  •  '  o 


,2Tr  /»27T  • 

/9m  n  =  /  sem(T7)  sen  (77)  di7/j  sen  (77^77 

.  }  .O'..,  “  • 

zw.  .:  //*27r 

Xm;n,  =  f  ce'jTj)  5^(77 )  *77/ ]q  sen*  ft)'*} 

-.6 

.27T  •*  //*27r 

Xm,n  -f  sem(77)cen  (77)  ^/J.  cen  (77)  d?7 


where  the  prime  signifies  the  derivative  of  the. function  with  respect  to 
its  argument.  The  above  integrals  can  easily  be  integrated  using  the 
series  expansions  of  Mathieu  functions  in  terms  of  trigonometric  functions 
(Ref.  5). 
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